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We study the effect of peculiar motion in weak gravitational lensing. We derive a fully relativistic 
formula for the cosmic shear and the convergence in a perturbed Friedmann Universe. We find a 
new contribution related to galaxies' peculiar velocities. This contribution does not affect cosmic 
shear in a measurable way, since it is of second order in the velocity. However, its effect on the 
convergence (and consequently on the magnification, which is a measurable quantity) is important, 
especially for redshifts z <1. As a consequence, peculiar motion modifies also the relation between 
the shear and the convergence. 
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I. INTRODUCTION 

Mapping the large-scale structure of the Universe is 
one of the most important current challenges for cosmol- 
ogy. Weak gravitational lensing represents a promising 
tool to achieve this goal, since it is directly sensitive to 
the distribution of matter in the Universe, independent of 
its nature (baryon, dark matter...). Gravitational lensing 
describes indeed the deflection of light rays from distant 
sources by the gravitational potential along the line of 
sight. It induces consequently a modification of the shape 
of the sources. This distortion of images contains infor- 
mation about the evolution of large-scale structure, i.e. 
about the geometry and dynamics of the Universe (see 
e.g. [1, 2] and references therein). Weak gravitational 
lensing can be divided in two parts: the shear, that dis- 
torts the shape of the source; and the convergence, that 
magnifies or demagnifies it. Both of these effects have 
already been measured. 

Cosmic shear is detected through the correlations it 
induces on the ellipticity of galaxies. It was measured 
for the first time in 2000, by four independent teams [3]. 
Since then, many other experiments have detected cosmic 
shear in random patches of the sky [4, 5]. In the next 
few years, weak lensing surveys, like CFHTLS [6], the 
Dark Energy Survey [7] and Pan-STARRS [8] will deliver 

accurate measurements (~ 1 % level) over large parts 
of the sky. In the further future, even more challenging 
experiments hke Euclid [9], LSST [10] and SNAP [11] arc 
planned. 

The other component of weak lensing, the convergence, 
can be detected through the modifications it induces on 
the galaxy (or quasar) number density over a given flux 
threshold [12]. The convergence (or more precisely the 
magnification) has already been robustly detected using 
quasar-galaxy correlations (see e.g. [13] and references 
therein). Moreover, recently [14] has highlighted the pos- 
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sibility to measure accurately the magnification autocor- 
relations with the Square Kilometer Array (SKA) [15]. 
A precision of 10% is expected at the beginning and 
^ 1 % later on. Hence the convergence provides an 
additional precise observational quantity, useful to con- 
strain cosmology. In order to make optimal use of this 
observational information, one needs to understand the 
underlying theory of weak lensing accurately. 

In this paper, we present a fully relativistic description 
of weak gravitational lensing. More precisely, we calcu- 
late the Jacobi map, that relates the surface of a galaxy 
to its angular image at the observer position, following 
the formalism presented in [16]. This map describes the 
distortion of a light beam by density perturbations along 
the geodesic between the source and the observer. The 
shear and the convergence are then extracted from this 
application. Our derivation differs from the standard one 
in two points. 

First, in the usual derivation the 4x4 Jacobi map is re- 
duced to a 2 X 2 matrix, called the magnification matrix, 
by assuming that the source and the image belong to the 
same two-dimensional subspace, normal to the photon 
direction and to the observer four-velocity. The conver- 
gence is then defined as the trace of the magnification 
matrix and the shear as the traceless part. However, the 
source four-velocity differs generally from the observer 
four- velocity. As a consequence the source and the image 
do not belong to the same two-dimensional subspace. In 
this paper we take into account this difference. We show 
that the Jacobi map can not be reduced to a 2 x 2 ma- 
trix but only to a 3 x 3 matrix, the third line describing 
the projection of the source plane into the observer plane. 
We establish that this generates an additional shear com- 
ponent. However since this effect is second order in the 
peculiar velocity of the source, it is not large enough to 
be observed. 

Secondly, in the usual derivation the shear and the con- 
vergence are expressed as functions of the source confor- 
mal time. However conformal time is not an observable 
quantity. Hence, in this paper we calculate the shear and 
the convergence as functions of the source redshift, which 
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is observable. We show that this modification generates 
new contributions to the convergence. Whereas most of 
those terms can be safely neglected with respect to the 
standard one, we establish that the contribution of the 
source peculiar velocity is important, especially for red- 
shifts z ^ 1. More particularly, for surveys in which the 
sources are situated at redshift 0.5, we expect a modi- 
fication of order 50 % relative to the standard results. 
For larger redshifts of the source the effect of peculiar 
motion decreases, whereas the standard term increases. 
However at redshift 1, we still expect the velocity term 
to be 1 % of the standard term, i.e. measurable by 
the SKA. Furthermore, we show that the transformation 
from conformal time to redshift does not affect the shear 
component. As a consequence, the relation between the 
shear and the convergence is modified by peculiar veloc- 
ities. 

The paper is organized as follows: in Sec. II we derive 
a general formula for the magnification matrix valid in 
(nearly) arbitrary geometries. In Sec. Ill we apply this 
formula to a perturbed Friedmann Universe. In Sec. IV, 
we investigate in detail the shear component of the mag- 
nification matrix. Finally, in Sec. V we calculate the con- 
vergence component and we determine its angular power 
spectrum. We investigate also the relation between cos- 
mic shear and convergence. 

Notation: We denote four-vectors with Greek indices, 
k". Three-dimensional vectors are denoted bold face k, 
or with Latin indices fc'. We use the metric signature 
(-,+,+,+)• 



II. MAGNIFICATION MATRIX 

We consider an inhomogeneous and anisotropic Uni- 
verse with geometry ds = (j^u dx^dx'^ . We are interested 
in the propagation of a light beam in this arbitrary space- 
time. We follow the derivation presented in [16]. We de- 
note by If the phase of the light beam. The wave vector 
is then given by — —Va^. We construct the deviation 
vector field iJx" connecting two neighboring rays. Since 
all the rays of the beam have the same phase, the devi- 
ation vector satisfies Sx^ka = 0. Furthermore it obeys 
the geodesic deviation equation [17] 

= R-.^sk^k'TSx' , (1) 

where A is an affine parameter along the geodesies, 
= fc^Va represents the covariant derivative along 
geodesies, and Rp^g is the Riemann tensor associated 
with the metric ds^. 

We consider the case of a light beam emitted by a 
galaxy at spacetime position S and received by an ob- 
server at O (see Fig. 1). 

We denote by vq the observer velocity and vs the 
source velocity. The photon energy measured at the 




FIG. 1: A light beam emitted by a galaxy at spacetime po- 
sition S and received by an observer at O. At the observer 
position, the plane normal to the observer four-velocity is in- 
dicated. 

source, respectively at the observer is 

cos = -ko.iXs)vs , (2) 
Wo = -ka{Xo)vo ■ (3) 

The solution of equation (1) is then given by [18] 

Sx'^iXs) = J'}ji\s)Se^{\o), (4) 

where (56i"(A) = —k/^VpSx^iX) . (5) 

UJo 

Here J'^(As) is the Jacobi map, that relates the deviation 
vector at the source Sx"{Xs) = Sxg to the angular vector 
at the observer 69" {Xq) = S9q. As we shall see, these 
two four-vectors belong to two different two-dimensional 
planes, orthogonal to the source peculiar velocity (re- 
spectively the observer peculiar velocity) and the photon 
direction at the source (respectively at the observer). 

Each ray can be parameterized by its affine parame- 
ter A and three other parameters that label the ray. 
This parameterization is not unique and one can there- 
fore make changes of the form 

y'^g^f), x = x + h{y^). (6) 

Under this reparameterization the connection vector Sx" 
transforms as [17] 

Si" = Sx" + k^Sh . (7) 

At the source, we can therefore choose a parameteri- 
zation such that SxgVsa = 0. Moreover, 6xgka{Xs) = 
induces Sx^nsa = 0, where ns is the photon direction 
at the source: ns = ^k{Xs) — vs- Hence Sxg lives in 
a two-dimensional subspace orthogonal to the source ve- 
locity and to the photon direction. It lies consequently 
in the plane of the galaxy. 

In addition to the choice of 5* and h in equa- 
tion (6), there is another degree of freedom corre- 
sponding to the scales of the affine parameter on the 
different rays. We can therefore require that at O, 
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voak°' = ''^Oa^jx' ~ ~^Oi for all rays. With this 
choice, we get (Sx"Va wo = 0, which gives (5^^QUOa = 0. 
Furthermore, we can easily show that S0"{\)ka{X) = 0, 
which implies SOqUq^ = 0, where uq is the photon 
direction at the observer: no — ^fc(Ao) ~ vq- Hence 
86q lives in a two-dimensional suSspace orthogonal to 
the observer velocity and to the photon direction. It lies 
consequently in the plane of the observer. 

Therefore the Jacobi map J"p{\s) in equation (4) de- 
termines how the surface of the galaxy represented by 
is deformed into the observer angular vector 56 q, 
during propagation in an arbitrary spacctimc. It maps 
hence the surface of the galaxy to its angular image at 
the observer. 

We now want to relate this map to the magnification 
matrix defined in [16]. At the observer we construct an 
orthonormal basis (Ef (Ao), (-^o), "o('^o)' ^o(-^o)) • 
We parallel transport this basis along the geodesic. The 
subspace defined by (i?"(A), £"2 (A)) is called the screen 
adapted to i'o(A) and fc"(A). 

We can write the deviation vector Sx'g and the angular 
vector 50q in this basis. Using the fact that Sx^ka = 0, 
we have 



Sx'^iX) = -a(A)£f(A)-6(W(A) 
-^o(A)[ng(A)+t;S(A)] 



(8) 



?o(A) 



= -a(A)Sf (A) - 6(A)S2«(A) - :^fc«(A) . 

u)o 

Moreover, 66Qka{Xo) = S9QVOa = implies 

66^ = -ei{\o)E^{Xo) - e2{Xo)E^{Xo) . (9) 
In this basis equation (4) becomes then 




JKXo,Xs) 



0^2 l(Ao), (10) 



where 



J}{Xo,Xs)= (11) 

( E,^{\s)J^{Xs)Ef{Xo) Ei^{Xs)J^{Xs)E^{Xo) 
E2o.{Xs)J^iXs)EfiXo) E2^{Xs)J1{Xs)Ei^{Xo) 
V k^{Xs)J1{Xs)E^AXo) k^iXs)J^{Xs)E^{Xo) 

J^(Ao, As) is proportional to the magnification matrix of 
lens theory, i.e. the gradient of the lens map. 

Indeed, the lens map relates angles at the observer do 
to angles at the source 63 



0o 



ds = eo + a. 



(12) 



The gradient of this map, A^j is called the magnification 
matrix 



A) = 6) + ^ 



For small angles, we can write 63 = A - Oo, and conse- 
quently 

$s = (Ao - Xs)ds = (Ao - As)^ • do • (14) 
This implies 

i(Ao,As) = (Ao-As)^. (15) 

Usually, the component ^o(As) of the connection vector 
is neglected, and equation (10) becomes 

|)(As)=(^D;.(Ao,A5)j-(^;)(Ao), (16) 

where Dj is the 2x2 submatrix of Jj. The shear and the 
convergence are then extracted from Z?'^- (see e.g. [16]). 

However, in general is not sufficient to relate ^5 
(or equivalently 9s) to 9o- The third line of the matrix 
J does indeed not vanish and therefore it plays a role in 
the deformation of the source plane. This is due to the 
fact that even if the connection vector at the source 
and the angular vector 60 at the observer arc both two- 
dimensional, they do not live in parallel planes. The 
connection vector |g belongs indeed to the plane normal 
to vs and ns, whereas the angular vector Go belongs to 
the plane normal to vo and no- If vs ^ vo, these two 
planes are different. 

Hence if we choose (^Ef{Xo), E^iXo)) as a basis 
on the observer plane, the parallel transported vectors 
{E"{Xs), E^iXs)) do not form a basis of the source 
plane. The connection vector will then have a com- 
ponent along A:", that we have to take into account 
when calculating the shear and convergence of galaxies. 
The difference between the two-dimensional plane of the 
source and the one of the observer comes mainly from 
the difference between the peculiar velocity of the source 
and the observer. There is an additional effect, due to 
parallel transport, which is however very small. Hence, 
if Vs = Vo, we can approximate the two planes as equal 
and we recover the 2x2 magnification matrix . 

In section III, we calculate explicitly the full matrix 
J in a perturbed Friedmann Universe (where in general 
Vs ^o) and we show that the third line of the matrix 
J induces an additional shear effect, which is however 
negligible since its amplitude is of second order in the 
source velocity. 



III. THE MAGNIFICATION MATRIX IN A 
PERTURBED FRIEDMANN UNIVERSE 

We now calculate explicitly the magnification matrix 
in a Friedmann Universe with scalar perturbations. In 
longitudinal (or Newtonian) gauge the metric is given by 



(13) g^.t.dx^dx'' = 0? [-(1 -I- 2*)d?7^ (1 - 'l^)^ijdx^ dx^\ 



(17) 
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For perfect fluids like ordinary matter, dark matter and 
radiation, the metric perturbations 5* and <i> are equal. 
Hence we assume in the sequel $ = ^f. We restrict 
ourselves to a spatially flat universe {K — 0), so that 

lij = Sij. 

Furthermore, since lightlike geodesies are not affected 
by conformal transformations, we perform the calculation 
in the metric 



g^^dx^dx'^ = -^gi^^dx^dx^ 



(18) 



We can then easily relate the magnification matrix in the 
two metrics by remembering that angles are not affected 
by conformal transformations, but distances scale with 
the conformal factor a. Hence 



J{>^o,>^s) = asJ(Ao, As) . 



(19) 



We calculate therefore the magnification matrix related 
to the metric 



g^vdx^dx^ 



{I + 2*)dry2 + (1 - 2<i)5ijdx'dx' , (20) 



and then we simply multiply the result by as to obtain 
the correct magnification matrix in a perturbed Fried- 
mann Universe. 

In [19], we have calculated the Jacobi map associated 
with the luminosity distance in a perturbed Friedmann 
Universe, which relates angles at the source 59^ to dis- 



tances at the observer 5xq 



(21) 



This situation is the mirror of the lensing situation con- 
sidered here (see equation (4) and Fig. 2). The two maps 
can therefore easily be related [17] by 

J"^(Ao, \s) = -— J^'"^/(As, Ao) . (22) 




Lensing 



Luminosity distance 



FIG. 2: On the left the "lensing" Jacobi map relates a surface 
at the source dAs to a solid angle at the observer dilo, and 
on the right the "luminosity distance" Jacobi map relates a 
surface at the observer dAo to a solid angle at the source 
dQ.s- The two maps are linked by equation (22). 



Using J 



(di) a 



from [19], we find 



/l-K-71 -72 0\ 
J){vo,Vs) = -asivo-Vs) -72 1 - k + 71 

\ w • El w ■ E2 / 
(23) 

with 



71 



d V"^ ^g)(^o !lld^d^i,(^ElEi~EiEi){24) 

r,s VO - VS 



72 = 2 / dri 
'vs 



vo {ji- ris)[rio - v) 



Vo - Vs 



d^d,-^E\Ei (25) 



dij V^W + 3ws 



vs 



Vo - VS 
mo 



ivc 



Vs) 



ns 



vo 



rf77*-2 / drj- — —ii (26) 



r,s VO - VS 



VS - Vo 



drjSI^ 



(27) 



Here V^^^" = (V^ — n^n^didj)"^ is the transverse Lapla- 
cian. Note that a similar derivation of 71 , 72 and k in the 
context of CMB lensing is presented in [20] . 

We have written here the magnification matrix as a 
function of conformal time 77. However, rj is not an ob- 
servable quantity. What we do measure is the redshift 
of the galaxy, which is also affected by perturbations, 
zs = zs + Szs- Now 



Jivs) = Jivi^s)) = J{zs) = J{zs) - -T^J{zs)5zs 

dzs 



Furthermore, 

jL ,Uzs) = ^ ( ^ 

dzs ^ + Zs xT-Csivo - Vs) 



and (see [19]) 

Szs = (l + zs) 

+ (vo - vs) • n 
Hence k becomes 



(28) 

— 1^ J + 1st order 
(29) 

rvo 

^'s - *o + 2 / dr]n ■ 

J vs 



(30) 



'"^ iv - Vs){vo - v) 
drj 



Vo - VS 



Vi* + 3*s 



2 / dr]^ — - ^ / d-q^ 

Ivs 'no - '7s [Vo - Vs) Jris 



1 - 



T-Lsivo - vs) 



(vs - Vo) • n 



+^s - *o + 2 / dr/* 



(31) 



and 71 , 72 and w are not affected. 
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The complete relativistie magnification matrix in a 
perturbed Friedmann Universe defined in equations (23)- 
(25), (27) and (31) contains therefore additional terms to 
the standard one. In the next section, we investigate 
the effect of the third line of the matrix on the shape 
of a galaxy, and more particularly on the shear. And in 
section V, we discuss in more detail the different con- 
tributions to the convergence k. We determine that the 
only term which can be relevant is the one involving the 
galaxy peculiar velocity. 



IV. THE SHEAR 

In this section, we study the effect of the third line of 
the magnification matrix on the shape of a galaxy. In 
order to simplify the calculation, we restrict ourselves to 
the peculiar velocity contribution. This means that we 
consider a homogeneous and isotropic Friedmann Uni- 
verse, but we allow for nonzero peculiar velocity of the 
source vg and of the observer vq- The magnification 
matrix becomes then 



Vo - r]s 



(32) 



1-K 







1 - K 



^ + ZS V (V5 - Vo)Ei (V5 - Vo)E2 , 



where 



K = 1 - 



Ti-sivo - Vs) 



(vs - vo) • n . (33) 



The usual shear components 71 and 72 vanish, but as 
we will see the third line generates an additional shear 
deformation. 

From equation (10) and (32), we find 



?0 = H 



1 — K 1 — K 



(34) 



where Vi = (vg — vq) • Ei and V2 = (vg — vq) • E2. At 
first order in vg and vq, equation (34) becomes 



Co - Vi^i - V2^2 = 



(35) 



that represents the equation of a two-dimensional plane. 

This reflects directly the fact that the galaxy, described 
by does not belong to the parallel transported ob- 
server plane: (E'f (Ag), E2{Xs)), but it rather belongs 
to the plane defined by equation (35). 

We now assume that the galaxy is a disc of radius r. 
Hence in spherical coordinate 

Ci = rsin^cos^ , ^2 = rsin^sin^ , = t"cos^ . (36) 

Combined with equation (35), this gives 



6(0) = ±r 



COS( 



sm ( 



(37) 



with t{(j}) 



v\ cos <p + V2 sm ( 



(38) 



The -I- sign is for ((> £ [0, 7r/2] U [tt, 37r/2] and the - sign 
for G [7r/2,7r] U [37r/2,27r]. 
The observer measures 



m 



92(0) = m=± ' 



smcj) , (39) 



where we have already removed the monopole contribu- 
tion 

We now show that (^i ((/>), ^2 (</>)) describes an ellipse, 
and we determine its semiaxis a and h. We consider, 
without loss of generality, the case where W2 = 0. This 
simply means that we align Ei on V5 — vo- We then 
determine a and b such that 



Using equation (39) for 9i{(j)), we find 



1-K 



and 



2\-l/2 



(40) 



(41) 



Hence we see that the third line of the magnification 
matrix J deforms a disc into an ellipse and consequently 
it induces an additional shear effect. Since this new con- 
tribution does not derive from a scalar potential, as the 
usual component, it can generate B-modes. Peculiar mo- 
tion constitutes therefore an intrinsic source of B-modes 
in cosmic shear. Moreover, contrary to the B-modes' con- 
tribution from source redshift clustering [21] that peaks 
at small scales, the velocity contribution is expected to 
peak at rather large scales, where peculiar velocity corre- 
lations are larger. However, this effect is second order in 
the velocity difference v\. It is therefore too small to be 
detected by current and future experiments and conse- 
quently to be responsible for the observed B-modes [22]. 
In the following we can therefore safely neglect it with 
respect to the standard shear components 71 and 72 and 
reduce the magnification matrix to the usually considered 
2x2 submatrix. 

Note that in this calculation we did not take into ac- 
count the effect of the potential, which induces an in- 
tegrated Sachs- Wolf term, J^^drjVif in equation (27). 
However, this term is much smaller than the peculiar ve- 
locity contribution (sec e.g. [19]) and can be neglected. 
Hence, the effect of the third line on the shape of the 
galaxy reduces to a purely kinematic effect, which can 
also be understood as length contraction of the galaxy in 
the velocity direction. 
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V. CONVERGENCE 
A. The velocity contribution to the convergence. 

Wc evaluate now the convergence k. The various terms 
in equation (31) are very similar to the one affecting 
the luminosity distance perturbations. In [19], wc have 
estimated all these terms and we found two dominant 
contributions: the lensing term and the source peculiar 
velocity term. Wc neglect hence the other terms here. 
Among them, one finds the effect of peculiar velocity of 
the lenses, which corresponds to the gradient of the po- 
tential in equations (27) and (30). Those terms which 
may be relevant in intermediate lensing (see e.g. [23] and 
references therein) are completely subdominant in weak 
lensing. Therefore in the following we restrict our calcu- 
lations to the two components 



where 



Jrt 



'VS 

Kv = I 1 - 



Vo - VS 
1 

'Hsivo -Vs) 



(42) 



(43) 



(vs - vo) • n . (44) 



The convergence can be measured through the modifica- 
tions it induces on the galaxy number density at a given 
fiux. 

Let us introduce the magnification 



det^ 



1 + 2k, if K,7i,72 < 1. (45) 



The magnification modifies the size of an observed source: 
dflo = f^dils, where dQs is the true angular size of the 
source and dflo is the solid angle measured by the ob- 
server, i.e. the size of the image. The lensing term is 
always positive and consequently it always magnifies the 
source. On the contrary, the velocity term can be ei- 
ther positive or negative and it therefore either magnifies 
or demagnifies the source. The sign of depends on the 



1 - 



sign of vs-n and on the sign of g{zs) = Hs(rio-vs) , 

To fix the ideas, let us consider a galaxy moving to- 
ward us. Then • n > 0, since n represents the photon 
direction, which points to the observer. Moreover, the 
sign of g{zs) depends on the redshift of the source zs- In 
a ACDM Universe with = 0.24 and flji = 0.76, we 
have g{z) < for 2; < 1.7 and g{z) > for 2; > 1.7. Hence 
for small z the surface is dcmagnified {p,^ = 1 + 2kv < 1) 
and for large z it is magnified (/Uv > !)• This is caused 
by the change in redshift induced by the source peculiar 
velocity. Indeed, for a fixed redshift, a source moving to- 
ward us is more distant (in conformal time for example) 
than a source with null peculiar velocity. This generates 
two opposite effects on the observed solid angle. On one 
hand, a distant galaxy is observed under a smaller solid 



angle. And on the other hand, since its conformal time is 
smaller, its scale factor is also smaller. Consequently, the 
image experiences more expansion when coming to us. At 
small redshift the first effect dominates, leading to a de- 
magnification of the source, whereas at large redshift the 
second effect dominates and the source is magnified. At 
z = 1.7, both effects compensate, leaving the size of the 
source unchanged. The situation is then simply reversed 
for a source moving away from us. 

We evaluate now the effect of magnification (or demag- 
nification) on the galaxy number density. We consider 
n{f)df unlcnscd galaxies per unit solid angle at a redshift 
zs and with a flux in the range [/, f + df]. The magnifi- 
cation modifies the flux measured by the observer, since 
it modifies the observed galaxy surface. It modifies also 
the solid angle of observation and hence the number of 
galaxy per unit of solid angle. These two effects combine 
to give a galaxy number overdensity [13] 

= "^^l7,f-^^ ^ 1 + 2(a - !)(«* + «v) . (46) 

Here a = -A^'(> fc)Ic/N{f,), where N{> /c) is the 
number of galaxies brighter than fc and fc is the fiux 
limit adopted. Hence a is an observable quantity [13, 14]. 

Recent measurements of the galaxy number overden- 
sity (5^ are reported in [13]. The challenge in those mea- 
surements is to eliminate intrinsic clustering of galaxies, 
which induces an overdensity 5'^ much larger than 5^. 
One possibility to separate these two effects is to cor- 
relate galaxy number overdensities at widely separated 
redshifts. One can then measure {5'^{zs)5'^{zs'))i where 
Zs is the redshift of the sources and zs' < zs is the red- 
shift of the lenses. Recently, the authors of [14] proposed 
to remove close pairs of galaxies at the same redshift 
from the signal in order to eliminate 5f . This allows 
then to measure {5'^{zs)5'^{zs')) , cither for zs 7^ or 
for Zs = zs', i-e. for galaxies situated at the same red- 
shift but in different directions. This method requires of 
course to know precisely the redshift of the galaxies. 

The velocity contribution Kv has only an effect on 
{5^{zs)S^{zs'))- The correlations between 5'^{zs') and 
V5 are indeed completely negligible and hence the source 
peculiar velocity does not affect {5'^{zs)5g\zs'))- In the 
following we study in detail the contribution of peculiar 
motion to {5i^g{zs)5^{zs')). 

The two components of the convergence and Kv 
(and consequently the galaxy number overdensity) are 
functions of redshift zs and direction of observation n. 
We can therefore determine the angular power spectrum 

{5>^g{zs, n)5,^(2s', n')) = ^^^^^^5, ^5')^^(n ' n') • 

(47) 

The coefficients C^{zs, zs>) contain two kinds of terms 
induced by {k^k^) and (kvKv)- The cross-term (kvK*) 

vanishes since Kq, contains only Fourier modes with a 
wave vector perpendicular to the line of sight (see 
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eq. (43)), whereas Kv selects modes with wave vector 
along the line of sight (eq. (44)). 
The velocity contribution is 



4(as - l){as' ~ l)((vsn)(vs/n')) • 
1 \ 1 



(48) 



1 - 



^5(^70 - ^75) 



• 1- 



T-Cs' ivo - VS' ) 



Here we neglect the peculiar velocity of the observer vo 
which gives rise only to a dipole. 

We use the Fourier transform convention 



v(x,r;) = ^ J d^kyr{k,rj)e^^- , 
v(k, '?) = y d^^^i^, ?7)e~*'"' . 
With the solution of the continuity equation [24] 



(49) 
(50) 

(51) 



where (5(k, 77) is the density contrast, D(a) is the growth 
fimction, and D'{a) its derivative with respect to 77, we 
find 



CJ{zs,zs') = 



16TT6j,{as - l){as' - 1) D' (zs)D' jzs') 



'Hsi'no - Vs) J \ 'Hs'ivo - VS') , 
rfo / dkkT\k)jmrio - mMiKm ~ m')) (52) 



Here 5h is the density contrast at horizon and T{k) is 
the transfer function defined through [24] 

*(k,a) = Avl,^(k)r(fc):^ , (53) 

where 'l'p(k) is the primordial power spectrum. 

We want to compare this contribution with the usual 
contribution coming from the potential 



Cf{zs,zs>) = 



D^{a = 1) 

D[zs)D{zs-){l + zs){l + zs')- I dkk^T\k) 



(54) 



VO 



Vs' 



dr)W{r),ris)(^jeik{r]o - v))+jeik{vo - v))) ' 
dv'W{v',Vs'){Mk{vo - v'))+JeiHvo - v'))) , 



where W{r],ris) = ^'^ (^qI!.^") '''' is the lensing kernel. 
We evaluate and Cf in a ACDM Universe with 
= 0.24, JIa = 0.76 and Ug = 1. Sh can 
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FIG. 3: The velocity contribution CJ as a function of I, for 
redshift (from top to bottom) z = 0.2, 0.5, 0.7, 1, 1.2 and 1.5. 
We see that CJ decreases with redshift in this redshift range. 



be measured from the Sachs- Wolf plateau [25]. We find 
Sh ^ Q - 10"^ [26]. We use the BBKS transfer func- 
tion [27] (see also [28] which proposes a similar fit). 

In Fig. 3-6, we plot CJ and Cf for different redshifts, 
but with zs — zs' ■ The amplitude of CJ and Cf depends 
on (a — 1)^, which varies with the redshift of the source, 
the flux threshold adopted, and the sky coverage of the 
experiment [14]. Since this term influences CJ and Cf 
in the same way we do not include it in our plot. Gen- 
erally, at small redshifts, (a — 1) is smaller than 1 and 
consequently the amplitude of both CJ and Cf is slightly 
reduced, whereas at large redshifts {a — 1) tends to be 
larger than 1 and to amplify CJ and Cf [14]. However, 
the general features of the curves and more importantly 
the ratio between CJ and Cf are not affected by (a — 1). 

Figures 3 and 4 show that CJ peaks at rather small £, 
between 30 and 150 depending on the redshift. The be- 
havior of cf is different, since it increases with i (Fig. 5 
and 6). Hence, peculiar motion of galaxies generates ad- 
ditional correlations that peak at relatively large angle 
6* ~ 70 — 360 arcmin. It is therefore important to have 
large sky surveys to detect this effect. 

The relative importance of CJ and Cf depends stron- 
gly on the redshift of the source. At small redshift, zs = 
0.2, the velocity contribution is about 5 ■ 10~^ and is 
hence larger than the lensing contribution which reaches 
10-6. At redshift zs = 0.5, CJ is about 50 % of Cf , 
whereas at redshift Z5 = 1 , it is about 1 % of Cf . Then 
at redshift zs = 1.5 and above, CJ becomes very small 

with respect to Cf : CJ ~ 10~^C*. This fast decrease of 
CJ is a consequence of the behavior of g^(z) plotted in 
Fig. 7. We see that at redshift zs = 1.7, CJ vanishes, due 
to the fact that = 0, as explained above. Finally, for 
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FIG. 4: The velocity contribution CJ as a function of I, for 
redshift (from top to bottom) z = 4, 3, 2.5 and 2. We see that 
CJ increases with redshift in this redshift range. 



FIG. 6: The potential contribution Cg as a function of I, for 
redshift (from top to bottom) 2 = 4, 3, 2.5 and 2. 




zs > 1.7, g{z) and consequently CJ start to increase as 
a function of zs (Fig. 4). However, the lensing term Cf 
increases faster than the velocity term and hence the ratio 
CJ /cf stays small, even at redshift 4 where it reaches a 
few 10-*. 

We have also calculated CJ{zs,zs') for zs ^ zs'. We 
found that it decreases very rapidly with the redshift dif- 
ference, and hence it gives a negligible contribution to 
correlations between different redshift bins. On the con- 
trary, cf (zs, zs') is not very sensitive to the redshift 



difference and induces therefore correlations at different 
redshifts. 

To summarize, we see from Fig. 3 and 4 that a lensing 
survey must satisfy three criteria in order to detect the 
velocity contribution CJ. First it has to cover a rather 
large part of the sky, since the velocity contribution peaks 
at angles between 70 arcmin and 6 degrees (depending on 
the redshift of the source). Ongoing and future surveys 
fulfill this requirement. For example, the CFHTLS Wide 
survey has already spanned 57 square degrees of the sky 
(over three independent fields) and allowed consequently 
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to measure the two-point shear statistics from 1 arcmin 
to 4 degrees. And in the near future, measurements will 
be extended up to 8 degrees [5] . Future surveys are even 
more ambitious. For cxampl(\ Euclid and the SKA plan 
to cover up to 20'000 square degrees. The requirement of 
a wide field survey is also of great importance for shear 
measurements. It is necessary in order to, on one hand, 
beat cosmic variance, and on the other hand probe accu- 
rately the linear regime of perturbations. 

Second, in order to measure precisely the galaxy num- 
ber overdensity, one needs to remove intrinsic clustering. 
One way to perform this, consists simply of removing 
from the signal close pairs of galaxies at the same red- 
shift, as proposed in [14]. This requires one to measure 
the redshift of the galaxies. The SKA offers here a great 
advantage, since it can measure with high precision the 
redshift of galaxies from 21cm emission line wavelength. 
Most other future lensing surveys will deliver photomet- 
ric redshift measurements of the observed galaxies. This 
is indeed also crucial for cosmic shear observation. First 
it allows, as for magnification measurements, to remove 
close pairs of galaxies, in order to reduce the systematic 
noise due to intrinsic alignment of galaxies in the same 
cluster [29]. Second, redshift information permits one to 
perform 3D cosmic shear analysis. 

Finally, Fig. 3 shows that the velocity contribution is 
mainly important at redshift 2: ~ 1. Hence an ideal sur- 
vey should cover redshifts from 0.2 to 1. This is well 
within the range of ongoing and future experiments. 

The accuracy with which the velocity contribution can 

actually be measured as well as the exact redshifts range 
where it is detectable depend of course of the precision 
of the survey. In [14], measurements of the magnification 
autocorrelations by the SKA are discussed in detail. The 

autocorrelations are calculated for redshifts z ~ 1. Since 
at redshift 1, our velocity contribution is of order 1% of 
the standard term, it should be measurable by the SKA 
if, as expected, an accuracy of 1% is reached. At larger 
redshifts, the velocity contribution decreases quickly and 
will hence not be detectable. However, the SKA should 
be able to measure the velocity contribution at smaller 
redshifts. Indeed, at small redshift z ~ 0.2, the velocity 
contribution is rather large, of the order 5 • 10~^. Hence 
between redshift 0.2 < 2; < 1, the velocity contribution 
should be measurable. Between redshift 0.5 <z< 1.5, a 
large number of galaxies 10^ deg~^) are expected to 
be detected by the SKA, hence good statistic should be 
achieved in this redshift range. Of course, the amplitude 
of the signal is weighted by (a — 1)^, which is not included 
in Fig. 3 and 4. At very small redshift, a — 1 can be small, 
and even negative. If (a — 1) ^ —0.5, the amplitude is 
for example reduced by 4. The value of a depends on the 
features of the survey: sky coverage, selection threshold, 
redshift of the galaxies, etc... Hence one can optimize 
this value to detect the velocity contribution. 



B. Nonlinecir contribution 

Until here, we have only considered peculiar veloci- 
ties from linear theory, where equation (51) holds. How- 
ever, on small scales, galaxies' peculiar velocities do 
not obey linear theory anymore. Typically, galaxies 

in a cluster can have peculiar velocities of the order 
vs ^ 1000 km/s [30]. Equation (48) becomes then 

(5;(zs, ^WAzs, n')) ^ {as - V)Hg\zs) ■ 10"^ . (55) 

The amplitude depends on the prefactor 4g^(z) plotted 
in Fig. 7, and on (a — 1)^. Hence the nonlinear contribu- 
tion is large up to redshift z ~ 1.4, where it is still of the 
order of 10~^, if we take a — 1 ~ 1. Then between z = 1.4 

and z = 2.1 it becomes very small. But for z ~ 2.1 the 
signal is again non-negligible and reaches even 10~^ at 
z = 4. However at those redshifts, we do not expect to 
find clusters of galaxies. 

Moreover, the nonlinear velocity contribution is large 
only on very small scales and it should become negligible 
when one removes close pairs of galaxies at the same 
redshift from the signal, as proposed in [14]. 

C. The reduced shesir 

In section IV we have established that peculiar motion 

contributes at second order to the shear and is therefore 
too small to be detected. However, since what is actually 
observed is not the shear, but the reduced shear [37] 

g = ~ 7 — K7 for K, 7 <C 1 , (56) 

i K 

it is necessary to investigate the effect of peculiar motion 
on g. 

The reduced shear correlations are given by 

{gizs, n)g(zs', n')) = (7(2^5, n)7(zs', n')) (57) 
-(7(2^5, n)7(zs', n')K(zs', n')) 
-(7(2^5, n)«;(zs,n)7(zs',n')) + higher order . 

The corrections (7(z5, n)7(zs', n')K5,(zs', n')) have been 
computed in [38]. They are large enough to impact on 
cosmological parameter estimation from future experi- 
ments. However, one can show that the velocity induced 
corrections (7(z5, n)7(zs', n')Kv(-2S') n')) vanish. Indeed 
7(zs,n) contains only Fourier modes with a wave vec- 
tor k_L perpendicular to the line of sight (see eq. (24) 
and (25)). On the other hand Kv(zs,n) selects modes 
with wave vector along the line of sight (eq. (44)). Hence 
correlations between and 7 vanish and we have 

(7(zs,n)7(zs',n')Kv(zs',n')) = (58) 
(7(^s,n)7(^S',n'))(Kv(^S',n')) = . (59) 

Therefore peculiar motion does not affect the reduced 
shear correlations at first order in Kv- 
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D. Relation between shear and convergence 



Similarly in the all sky calculation, we have 



In the usual description of weak lensing, the shear and 
the convergence satisfy a functional relation. They are 
indeed both related to second order derivatives of the 
potential along the line of sight. In the flat sky approxi- 
mation, we have 

K*(r/s,/3) = Afs^{r]s,/3) , (60) 

7i(r?s,/3) = (V/3,V/3, -V;3,V/3,)*(r?s,/9) , (61) 

72(r/s,/3) = 2V/3,V;3,*(r?s,/3) , (62) 

where 

rvo 

^(j?s,/3)=/ dvW{v,7js)^iv,(3) . (63) 

Here /3 = (/3i,/32) is a two-dimensional angular vector 
describing the position of the source in the flat sky ap- 
proximation, and V/3 is the associated two-dimensional 
gradient. Equation (60) can be solved to express as a 
function of k\s,. Inserting the result into equations (61) 
and (62) leads to the following relation between and 
7 (see e.g. [31]) 

7(/3) = 71 (/3) + »72(/3) = d^P'D{(3 - (3')n^(P') , 



where 



|/3|4 



(64) 
(65) 



The two-dinicmsional Fourier transform of the shear and 
the convergence are therefore related through 7(q) = 
e^'^K*(q), where <^ is the polar angle of q. Consequently, 
the two-dimensional power spectrum of the convergence, 
Pk^ and of the shear, are equal. 

A similar relation holds in the all sky calculation, 
where the angular power spectra satisfy [32] 



(€ + 2)(£-l) ^ 



(66) 



These relations are modified by the velocity contri- 
bution. Galaxies' peculiar velocities generate indeed a 
new contribution to the convergence, /tv that is not re- 
lated to the potential ^ . On the other hand, it does 
not change the shear at flrst order. Hence, in the flat 
sky approximation the measured shear 7 and convergence 
K:obs = + Kv obey 

7(/3) = cff3'D{^ - ^'){Kobs{0') - Kv(/3')) • (67) 

Since the correlations between and Kq, vanish, the 
power spectra satisfy 



(^ + 2)(f-l) 



(69) 

4(a - 1)2 ' ^^^^ 



(68) 



where we have used equation (46) for the second equality 
sign. 

Consequently, if one measures both the shear CJ and 
the magnification C'^ as functions of the redshift, equa- 
tions (69) allows one to extract the peculiar velocity con- 
tribution CJ . This provides a new way to measure di- 
rectly peculiar velocities of galaxies. Future surveys, like 
Euclid and the SKA should be able to deliver measure- 
ments of both the shear and the magnification in the 
redshift range 0.2 < z < 1. Since at those redshifts 
the difference between the velocity contribution and the 
potential contribution is at least of the percent level, it 
seems feasible to extract peculiar velocities from the mea- 
surements of and . 

A direct measure of the peculiar motion of galaxies 
would be extremely interesting, since it allows one to de- 
termine the underlying matter power spectrum. More- 
over, since the velocity contribution CJ peaks at rather 
large angles, it provides a measurement of the velocity 
field well inside the linear regime, where the relation be- 
tween peculiar velocity and matter overdensity is sim- 
ply given by equation (51). A more careful analysis 
is of course necessary to determine if this method can 
compete with actual or future observations of peculiar 
velocities, from galaxy surveys (see e.g. [33] and refer- 
ences therein), supernova surveys [34] or from the kinetic 
Sunyaev-Zel'dovich effect in the Cosmic Microwave Back- 
ground (see e.g. [35, 36] and references therein). 



VI. CONCLUSIONS 

In this work we have studied the effect of peculiar ve- 
locity on weak gravitational lensing. We have derived a 
general formula for the shear and the convergence, that 
takes into account all relativistic effects of linear pertur- 
bation theory. We have identified a new important con- 
tribution generated by the peculiar motion of galaxies. 

We have shown that the shear component is affected 
only at second order by peculiar velocity. Consequently 
this contribution does not aff'ect cosmic shear in a mea- 
surable way. However, we have found that the effect of 
peculiar velocity on the convergence is important and 
cannot be neglected, especially for redshift z ^ 1. At 
small redshifts z ^ 0.4, the peculiar velocity contribu- 
tion is even larger than the lensing contribution. Hence 
measurements of the convergence (or more particularly 
of the magnification) at small and intermediate redshifts 
are affected by the peculiar motion of galaxies. In the 
redshift range 0.2 < z < 1 we expect an effect large 
enough to be detected by future galaxy surveys, like, for 
example, the SKA and Euclid. 
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We have also shown that pecuhar motion modifies the 
relation between cosmic shear and convergence. One way 
to measure the effect of peculiar velocity would hence 
consist in measuring both the shear and the convergence 
at the same rcdshifts, and in comparing them to extract 
the velocity contribution. This provides hence a new way 
of measuring the velocity field directly. 

Finally, we have found that the redshift dependence of 
the velocity term is mainly given by the prcfactor g'^{z), 
plotted in Fig. 7. Since this function depends directly 
on cosmological parameters, it would be very interesting 
to find a way to extract it from the data, in order to 



constrain cosmological models. 
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